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Abstract 

The paper describes experiments with the SND detector at VEPP- 
2M collider, carried out during the period from October 1995 until June 
1997. The total integrated luminosity of 6.4 pb" 1 was collected in the 
energy range 2E — 0.4 -r 1.4 GeV (MHAD97 experiment), corresponding 
to 4 • 10 5 /i + /i~-pairs produced. Preliminary results of the 1996 <^>-meson 
experiment (FI96) are presented. The total number of ^-mesons produced 
is 4 • 10 6 . New data on rare decays of <fi and r/(550) mesons, in particular 
B((p -> 777) = (1-30 ± 0.06 ± 0.07) %, 

B(<j> 7r°7r°7) = (1.1 ± 0.2) ■ 10" 4 , (M^o > 800MeV), 
Bfo-»/°7) =(4.7±1.0)-10" 4 , 
B{<j}-^rfK° 1 ) = (1.3 ±0.5) ■ 10" 4 , 
B^^rj'-y) < 1.7- 10" 4 , 
B{(j>^2-K ) <6-10" 4 
were obtained. 
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1 The MHAD97 experiment 



In the end of 1996 the experiments, started in 1995 with SND detector ||, 
were continued. In the period from October until November the energy range 
2E from 600 down to 370 MeV was scanned. The total integrated luminosity 
of 80 n& -1 was collected, which corresponds to approximately 3 • 10 4 /i + //~ and 
10 4 7r + 7r - events. The goal of the experiment was to collect data needed to 
improve accuracy of calculations of the contribution of this energy region into 
anomalous magnetic moment of muon. The data were taken in parallel with the 
CMD-2 B detector. The energy dependence of VEPP-2M average luminosity 
during this experiment is shown in Fig. [j]. 

In the beginning of January 1997 the VEPP-2M energy was set to approxi- 
mately 2E = M,j, and data were collected corresponding to integrated luminosity 
of about 200 nb^ 1 or 0.7 • 10 6 ^-mesons produced. 

From the end of January until June 1997, in parallel with the CMD- 2 de- 
tector, the MHAD97 experiment in the energy range 2E from 980 to 1380 MeV 
was carried out. Two scans of this energy region were performed: one - up- 
wards, and another - downwards, with a step of A(2E) = 10 MeV. The 
total integrated luminosity collected was 6.3 pb" 1 at an average luminosity of 
1.3- 10 30 cto~ 2 s _1 (Fig. Q), and the total number of recorded events of 10 8 , from 
which about 3.6 • 10 were pairs, and 10 7r + 7r~- pairs. In Fig. | the 

weekly schedule of data taking is shown. The mean event recording rate was 
about ~ 23 Hz and average live time ~ 72 %. The raw events were stored on 
thirty 8 mm 4 GB tapes. 

The main goal of this experiment was a thorough measurement of dif- 
ferent hadronic production cross sections, including e + e~ — * 2tt 1 3-7T, 47r, W7r, 
KK, KKtt, .... New measurements can facilitate more precise tests of differ- 
ent theoretical models (VDM, CVC, ...), check if there are sizable contributions 
from radial excitations p',u>',(j)', and determine the contribution of this energy 
region into muon anomalous magnetic moment. 



2 Detector performance 

The SND detector was described in detail elsewhere in Jl], [|, ^]. Simultane- 
ously with data taking the work on improvement of detector parameters was 
performed. One of the most important parameters is the energy resolution of 
the calorimeter for electrons and photons. Primary calibration of the calorime- 
ter was done using cosmic muons, and final calibration - using Bhabha events. 
Recently, progress in understanding of factors limiting the detector resolution 
(Fig.||) was achieved. It turned out, that nonuniformity of the light collec- 
tion efficiency in the Nal(Tl) crystals of the inner calorimeter layer contributed 
significantly into calorimeter response. When this effect was included into 
simulation program, the disagreement between measured (<je/E 7 — 5.5% at 
E 1 = 500 MeV) and simulated resolution (ag/E^ = 4.2%) was significantly 
reduced. The next step in improvement of calorimeter response was implemen- 
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tation of an absolute calibration procedure based on Bhabha scattering events. 
With such a calibration the energy resolution for electrons and photons was 
improved by approximately 10%. Such a modest improvement with respect to 
" cosmic" calibration [^| shows relatively high precision of the cosmic calibration 
procedure. 

Another important parameter of the calorimeter is an invariant mass reso- 
lution a m /m for particles decaying into photons, which depends on a particle 
energy and mass. Table [l] shows values of a m /m for ir°, 77, Ks uj- mesons, 
extracted from FI96 experimental data || . 



Table 1: Invariant mass resolution (er m /m) for particles decaying into photons. 



Particle 


7r° 


V 


V 


K S 




Energy, MeV 


519 


657 


657 


510 


801 


Process 


<t> —* 7T°7 


4, -» rn 


4> -* irr 


<j> -> K S K L 




Decay mode 


7T° — * 77 


ri -> 77 


r) -> 3tt° 


K s -> 2tt° 


CJ — + 7T°7 


a m /m, MC, % 


8.1 


2.7 


2.9 


4.4 


1.8 


CT m /m, exp., % 


10 


2.9 


3.5 


5.2 


2.4 



Systematic discrepancy between measured and simulated resolutions, seen 
in the Table [l[ is still under investigation. For example, in Fig. [| the u>- meson 
invariant mass spectrum in the decay channel u> — » vt°7 is shown. It was obtained 
during the study of the process e + e~ — > lutt° — * 7r°7r°7. 

Before the beginning of MHAD97 experiment the performance of the drift 
chamber system was improved, using for calibration copiously produced collincar 
e + e~- pairs from Bhabha scattering process. The spatial resolution in transverse 
direction, determined by the accuracy of drift time measurements was equal 
to 180 microns (a) on average; the resolution in longitudinal direction was 
3.2 mm (a) for charge division, and 1.5 mm for cathode strip readout. In 
Fig. ^|, ^| the experimental and simulated distributions over Atp, Ai9 angles 
in the Bhabha scattering events are shown. The widths of the distributions 
are determined not only by spatial resolution of the drift chambers, but also 
by higher order QED corrections to Bhabha scattering. The intrinsic angular- 
resolution of the drift chambers is a,p = 0.7° in azimuth and a$ — 2.2° in 
polar plane. The performance of drift chambers can also be described in terms 
of impact parameter AR - measured distance between the track of a charged 
particle and beam axis, projected onto a plane, perpendicular to the beam 
(Fig. 0). In the infinite spatial resolution limit the average AR value must be 
equal to the transverse size of the beam a x = 10 microns, a r = 200 microns. 
For the Bhabha scattering events <7a.r — 0.4 mm was observed. 

The luminosity was monitored during the experiment using the processes of 
double bremsstrahlung (by an external small-angle monitor), Bhabha scattering, 
and 27annihilation. Listed below are the selection criteria for events of the two 
latter processes, used for luminosity monitoring. 

For the e + e~ — * e + e~ process: number of particles < 4; number of charged 
particles > 2; distance between tracks and beam collision point AR < 1 cm, 
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AZ < 10 cm; acollinearity angle A(p < 10°, A-d < 25°; polar angles of both 
tracks are within 27 < i) < 153° interval; energy deposition relative to the beam 
energy 0.6 < AE/E < 1.5 for both particles; two-charged-particles trigger must 
be produced in the event; 

For the e + e~ — > 77 process: number of particles < 4; number of neutral 
particles > 2; acollinearity angle between 7- quanta A<p < 10°, AD < 25°; polar 
angles of both photons are within 27 < < 153° interval; energy deposition 
relative to the beam energy 0.6 < AE/E < 1.5 for both photons; Two-neutral- 
particles trigger must be produced in the event. 

The luminosity is determined by the following expression for each process: 
Ni = Li<j Q i , where i is the fixed beam energy point number; iVj- number of 
e + e~ — * e + e~ or e + e~ — > 77 events; a i- detection cross section, calculated 
with an accuracy up to the third order in a. The ratio of luminosities i 77 /i ee 
for the third scan is shown in Fig. ||. In part of FI96 experiment Q this ratio 
was biased from its ideal value by several percent due to low gain in the inner 
drift chamber. 



3 Physical results of FI96 experiment 

The FI96 experiment itself was described in our previous preprint Q. In total, 
six successive scans were performed (F79601-^F79606) with the total integrated 
luminosity of ~ 4 pb^ 1 , corresponding to approximately 6.5 million (f>- mesons 
produced. At present, more than a half of the total recorded data are available 
for analysis (scans FY9602 F/9604). 

3.1 General description of the SND data analysis 

The SND data processing procedure consists of several successive steps || . 

1. Reconstruction of events, stored on primary tapes. For each event a list 
of particles with their parameters, including energies, angles, etc., is built. 
Reconstructed events, together with the parameter lists, are written to 
secondary tapes. 

2. Scanning of secondary tapes and creation of third-generation tapes, con- 
taining certain event classes, e.g. events containing 6-7 photons, or events 
with two charged particles and three photons, etc. 

3. Events of the classes, mentioned above are downloaded onto hard disks 
for further processing using well known application packages PAW [Q and 
MINUIT g, together with packages, developed in BINP: COCHA data 
management system [§, GIST histogram code fiojj, U NIMOD Monte Carlo 
simulation code jllj , ART tape archiving system [|12| , and other programs. 

The GIST histogram package calculates parameters, specific for the SND. 
Let us list some of them, the most widely used in data processing: 
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Etot/2E — total energy deposition in the calorimeter, divided by the center 
of mass energy 2E; 

E np /2E — sum energy deposition in the calorimeter by neutral particles, 
measured in the units of the center of mass energy 2E; 

Pt t/2E — absolute value of the sum of all particles momenta in an event 
normalized by a center of mass energy; 

7V 7 — number of photons, found in an event; 

N cp — number of charged particles; 

di — polar angle of the z-th particle (particles are sorted in the following 
way: charged particles first, then neutral particles, with descending order in 
energy within each group); 

$min — minimal spatial angle between particles and beam axis; 

(fi — azimuth angle of z-th particle; 

Ri — distance between the i-th particle track and beam axis in R— tp- plane; 
Zi — Z - coordinate of the i-th particle track intersection with beam axis in 
R — Z- plane; 

Xe — parameter, characterizing quantitatively the degree to which the 
energy- momentum conservation law is held in an event, analogous to x 2 m 
statistics; 

Xm — parameter, characterizing the degree of likelihood of assumption, that 
there are intermediate tt° or rj- mesons in an event; 

X 2 — parameter characterizing quality of photons in an event ]l3[ ] , ]l4| . 

The total number of parameters available under GIST package is more than 
2000. The parameters from this set could be further combined to produce 
complex constraints and logical formulas for event selection. 

During study of processes with substantial statistics the distributions of 
selected events over beam energy were fitted using the following expression: 

Nt = L{Ef) ■ {a{Ei) ■ 8{Ei) ■ 5 beam (E t ) ■ e{E t ) + a B {Bk)), (1) 

where is a number of events in the i-th energy point L — integrated 
luminosity in this energy point; a — theoretical cross section of the process 
under study; 6 — factor, taking into account radiative corrections, calculated for 
each energy point by convolution of theoretical energy dependence of the cross 
section with the probability for a photon with a certain energy to be emitted |15| ; 
Sbeam — factor, accounting for beam energy spread, also obtained by convolution 
of radiative corrected cross section with a Gaussian energy distribution of the 
beam particles; e — detection efficiency for a process under study for actual 
event selection criteria, calculated using Monte Carlo simulation; gb — total 
visible cross section for all background processes, calculated in the same way as 
for the process under study. 

Also were taken into account beam energy corrections in individual scans. 

The dependence of theoretical cross section on the beam energy for the 
processes proceeding via intermediate vector resonances V = p, u>, 4> was a P~ 
proximated as follows (see also pc|): 
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<E) = \ Y: V^- Dv{s) I 2 (2) 

Here a D v is a peak cross section; mv, Ty(s) are mass and width of the 
resonance respectively, 

D v (s) =m\-s- iy/sT v (s); s = AE 2 . 

The formulae describing the energy dependence of the resonance widths 
(IV(s)) and corresponding references are cited in |po| . 

To describe the interference between resonant and nonresonant processes the 
following formula was used: 

a(E) = a (l + a'(2E - m<t> )) ■ |1 - Z^^| 2 , (3) 

where a' is a derivative of the nonresonant cross section over energy; Z is a 
complex amplitude of the resonant process. 



3.2 Measurements of 0- meson parameters in the process 

— > K$Kl — > neutral particles 

The main goal of the study was to obtain the values of m^, T^, a a , and B((f> — > 
KsKl) with lowest possible statistical error in order to use them in the analyses 
of other processes with smaller statistics and for evaluation of systematic errors. 
Given the integrated luminosity of 3.4 ptr 1 and detection efficiency of e ~ 
15%, the number of detected 4> — > KsKl — > neutral particles events could be 
estimated to be equal to ~ 10 5 , which leads to a statistical error in a resonance 
mass: ~ T^^y 7 ]^ ~ O-OlMel^. Since the FI96 experiment consisted of six 
successive scans of the 4>- region, the study of the KsKl production makes 
possible to check the energy calibration of the collider. Usually the energy 
scale of the collider depends on many factors. To a relatively good accuracy 
of better than 20 keV the beam energy can be expressed in the following way: 
E(B, T) — k-B + r- T + C, where B and T are the magnetic field strength and 
the collider temperature, k, r, C are the parameters to be determined during 
collider beam energy calibration. The correction accounting for variations of 
the collider beam circulation frequency AE/E — —6 • A/// was also included. 

Events were selected with iV 7 > 4, 0.4 < E to t/2E < 1.2, with detected 
Ks — > 2ir° decay, and measured m,K s within 400 600Mey interval (Fig. ^). 
Detection efficiency for the <fr — > KsKl decay is equal to e = 14.8%. For 
the background — > 777 process, imitating the process under study due to 
misidentification of photons, it is equal to e\ = 1.1%. Experimental data were 
approximated as a sum of the process under study, resonant background from 
4> — > 777, and nonresonant background. The integrated luminosity was measured 
using e + e~ — > 77 process. The resonance excitation curve is shown in Fig. [To] 
and the results of processing of three scans are listed in the Table |[ 

Main conclusions from the Table ||: 
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Table 2: <f>- meson parameters from the e + e — > (f> — > K$Kl process. 





FI9602 


FI9603 


FI9604 


PDG,1996 


m f MeV 

cr , nb 
B(<j>^ K S K L ),% 
Number of events 

X 2 /N D 


1019.48 ± 0.03 
4.02 ±0.10 
1314 ±17 
30.1 ±0.4 
51247 
10/11 


1019.30 ±0.03 
4.59 ± 0.08 
1322 ± 13 
30.3 ±0.3 
78335 
43/12 


1019.20 ± 0.03 
4.10 ±0.08 
1405 ± 14 
32.2 ±0.4 
49871 
12/13 


1019.413 ± 0.008 
4.43 ± 0.05 
1488 ± 22 
34.1 ±0.5 



1. There exists a discrepancy between (j)- meson mass estimations in the 
second and fourth scans: Am = 0.28 ± O.OAMeV, which cause is yet 
unclear. 

2. In both 2-nd and 4-th scans x 2 /^d parameter is small, which is an argu- 
ment in favor of a good beam energy stability during these scans, at least 
in the vicinity of the (f>- meson peak. 

3. A beam energy leap presumably occurred in the FI9603 scan. This may 
explain larger <f>- meson width and high % 2 /Nd ratio in this scan. 

4. Taking into account the former, only FI9602 and FI9604 scans were used 
for (j)- meson width determination. The combined result is = 4.06 ± 
0.08AfeV '. This is by 3 standard deviations smaller than its table value. 
The corresponding values of decay parameters are the following: 

cr = (1360±ll±45)rab 

B{<j) -> K S K L ) = (31.0 ± 0.3 ± 1.0)%. 

We hope, that processing of the next three scans will clarify situation with 
systematic errors due to collider energy stability and detector performance. 

3.3 Radiative decays 
3.3.1 The cf) — > 777 decay 

The (/> — > 777 decay is a classic magnetic dipole transition from <fi- into 77- meson. 
In terms of nonrelativistic quark model |l9| this decay can be described as a 
quark spin flip inside 4>- meson: Si — > Sq + 7. The theoretical estimation 
and experimental value agree within about 10%. Up to now, more than 10 
measurements of 4> — > 777 branching ratio were done and its current table value 
is equal to 1.26 ±0.06%. In this paper preliminary results of new measurements 
of this parameter in three 77- meson decay modes are presented: 77 — > 27 (i?(?7 — * 
2 7 ) = 39%); 77 -> tt+tt-tt (B(r] tt+tt-tt ) = 23%) 77 3tt° (B(r) -> 3tt°) = 
32%). 

The cf> — > 777 — > 37 channel 
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The following event selection criteria were used: iV 7 = 3, 4; 0.8 < E to t/2E < 
1.1; P±/2E < 0.175; d min > 24°; E imm > 0.2E; 550 < m 77 < SOOMeV; 
X% < 30. Events from two <fi- meson scans were processed. Simulated events of 
the 4> — * ?77; 4> ~~ > KsKl', 4> —> ^I'i were processed in the same way in order to 
evaluate their detection efficiencies, which turned out to be: 22%, 0.01%, and 



8% respectively. In Fig. 11 the -q- meson invariant mass distribution is shown. 

The cross section energy dependence was approximated according to formu- 
lae from 0, H with a resonant background from the </> — > KsKl,tt°^ processes 
and a nonresonant background taken into account. The <p- resonance excitation 
curve in the </> — > 777 channel is shown in Fig. [j^. The values of fit parameters 
are listed in Table ||. 



Table 3: Main cf>- meson parameters measured using the e + e — * <j> — ► 777 — ► 37 
process. 



Experiment 


FI9602 


FI9604 


PDG,1996 


m,/,, MeV 
V^MeV 
cr (e+e _ — > tfc — * 177), nb 

177), % 
Number of events 

X 2 /N D 


1019.48 

4.02 
58.9±1.9 
1.34±0.05 
4256 
16/10 


1019.2 
4.1 
59.6±2.2 
1.40±0.05 
4040 
20/12 


1019.413 ±0.008 
4.43 ± 0.05 
53.2±2.5 
1.26 ± 0.06 



Combining data from both scans one can obtain the following branching 
ratio: 

B{(f> -► 777) = (1.37 ± 0.04 ± 0.08)%, 

where the first error is a statistical one and the second is systematic. 

The (j> — > 777 — > 37r°7 channel 

Due to high efficiency of the SND calorimeter to multiphoton events, full 
reconstruction of </>—!■ 777 —> 37r°7 — > 6, 77 turned out to be feasible. The char- 
acteristic feature of this channel is a peak at 77- meson mass in the recoil spec- 
trum of the most energetic photon in an event (Fig. |l3|). Events were selected 
according to following criteria: N 7 = 6, 7; E tot /2E = 0.8-=- 1.2; P to t/2E < 0.12; 
flmin > 27°, Xe < 25, x 7 < 20. For the energy dependence of the cross section 
the standard approximation was used (§). The energy dependence of the de- 
tection efficiency was neglected. The background of ~ 2 % due to the process 
4> — > KsKl was subtracted using the events outside of the main peak in the 
Fig. H 

In order to estimate SND systematic errors for multiphoton events with full 
reconstruction of all photons, an independent analysis of events with iV 7 = 7 
was carried out. In these events the presence of three 7r°- mesons was required. 
The detection efficiency in this case was lower (Table ^), but no background 
from other (f>- meson decays was observed at a present level of experimental 
statistics. 

It is clearly seen from the Table ^ that the cross sections and branching 
ratios of the process <f> — > 777 — > 37r°7 for the events with 7V 7 = 6, 7 and _/V 7 = 7 
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Table 4: Results of the study of 4> — > 777 — > 7r 7r°7r°7 decay for three independent 
scans of the 0- resonance region. 





FI9602 


F19603 


FI9604 


Number of events 6, 77 


1167 


1748 


1036 


N-j 1 /N§ 1 (experiment) 


0.60 ±0.04 


0.62 ± 0.04 


0.62 ±0.05 


Nfj/N^j (simulation) 


0.60 ±0.02 


0.60 ±0.02 


0.59 ±0.02 


Number of fully 77 


360 


564 


315 


reconstructed events 








Efficiency 6, 77 


8.0 ±0.1 


7.9 ±0.1 


7.7±0.1 


Efficiency 77 


2.72 ±0.07 


2.64 ±0.07 


2.57 ±0.07 


B(tj> -> 777) 6, 77, % 


1.24 ±0.04 


1.23 ±0.03 


1.31 ±0.04 


£(</> -» J77) 77, % 


1.14 ±0.06 


1.17 ±0.05 


1.20 ±0.07 



differ by approximately 8% (2a). This systematic error can be explained by 
misidentification of photons in 7 gamma events. 

The data for iV 7 =6,7 from all three scans were combined and the following 
final result was obtained: 

a (4> -► 7/7 3vr°7) = (54.6 ± 1.3 ± 3.5)n6 

B(<j> -> 777) = (1.25 ± 0.03 ± 0.08)%. 

The (j> — > 717 — * 7r + 7r _ 7r°7 channel 

Since this decay channel in addition to three photons contains two charged 
particles in its final state, the performance of the drift chambers can contribute 
into systematic error of the result. Comparison of this result with those for pure 
neutral channels can be a good probe of the tracking system performance. 

Events were selected according to following criteria: N cp — 2, N 7 > 3, spa- 
tial angle between charged particles 0^2 < 150°. Events must satisfy energy 
and momentum conservation laws with additional requirement, that two photon 
pairs must have an invariant masses close to that of tt°. This requirement was 
imposed within a kinematic fitting procedure Q. This procedure calculated 
the following parameters: the minimum of logarithmic likelihood function J^-y, 
— the r\- meson mass estimate, Dl — difference between minimums of 
logarithmic likelihood functions for two best combinations of photons, produc- 
ing required 7r°-s, — difference in likelihood functions minimums in 777 
and 7r + 7T~7r° hypothesis. In Fig. [l4] experimental and simulated distributions 
over reconstructed mass of the 77-meson are shown. After addition of two more 
requirements: L V1 < 15, 500 < m n < 600MeV, the detection efficiency was 
estimated to be 30.4 ± 0.9%, or, taking into account B(r) — ► 7r + 7r~7r°) = 0.23, 
was equal to £(777 — > 7r + 7r~7r°7) = (7.0 ± 0.2)% 

Processing of the FI9602 scan gave the following results on peak cross section 
and decay branching ratio: 

cr (e + e~ -> (\> -> 777) = 42.7 ± 3.1 ± AAnb 

b(<(> -> m) = (°- 98 ± °-° 7 ± 0.10)%. 

Only background subtraction error and statistical error of simulation were 
included into the systematic error. 
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The data were also processed under more stringent selection criteria: N cp = 
2, N 7 > 3, ai, 2 < 130°, R h2 < 0.5 , L I)7 < 10, D L > 2, D L3n < -0.5, 500 < 
m v < 600 MeV . The detection efficiency here is smaller, e(?77 — > 7r + 7r _ 7r°7) = 
5.4±0.2%, but the background was additionally reduced by a factor of 3. Thus, 
the background subtraction error here was significantly lower. 

o- (e+e- -> cj) -> 777) = (48.8 ± 3.8 ± 3.8)n& 

B(cf> -> J77) = (1.04 ± 0.08 ± 0.08)%. 

Since the systematic error of the latter result is estimated to be smaller, we 
consider the last result as a conclusive one for this channel. 

Conclusions on the (j> — > 777 channel 

Results obtained for these three different channels do not agree well. Ad- 
ditional study of systematic errors is necessary. Taking into account these dif- 
ferences and common systematic errors, the following combined results can be 
obtained: 

<7 (e + e- 777) = (56.8 ± 2.6 ± 3.1)n6 
B(<f> -> 777) = (1.30 ± 0.06 ± 0.07)%. 

3.3.2 The 4> -> tt°7 decay 

The <j) ~^ Wl decay was studied on the basis of one scan of 4> -meson region, 
corresponding to a total integrated luminosity of 0.4p& -1 i.e. 10% of experimental 
statistics. The preliminary result was already published in || and is equal to 
B{(j> 7r° 7 ) = (0.10 ± 0.02)%. 

3.4 Rare radiative decays 
3.4.1 Search for the cf) — > 7/7 decay 

Radiative decay <f> — > 7/7 is a good probe of an internal structure of the rj'- 
meson @ (jj, 

As it was mentioned in JT^ ], the (j> — ► 7/7 decay branching ratio strongly 
depends on the gluonium and strange quarks contents in 77'- meson. In an 
assumption of an ordinary quark structure of 77', predicted branching ratio is 
B{cj> -> 7/7) = 7- 10~ 5 . 

Two attempts to measure this decay were done earlier. The first upper limit 
was obtained with Neutral Detector ^0|: B{<p — > 7/7) < 4.1 • 10~ 4 ; and recently 
was published a preliminary result from CMD-2 detector pi] : B(<p — > rf^f) < 
2.4 • 10~ 4 . 

In this work the process 

<fi — > 77^ — * 7r + 7r~?77 — * 7T + 7r _ 7r + 7r _ 7r°7 — * 7r + 7r~7r + 7r~37 (4) 
was studied. 

Among all decay modes of 77'- meson this one has a relatively high prob- 
ability: -B(t7' — * ir + ir~ri) ■ B(rj — ► 7T + 7r~7r°) ~ 10%, well defined final state, 
practically no combinatorial background in neutral particles, and it also con- 
tains an 77 meson in the decay chain. 
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Event selection for the search of the process (Q) was done using the following 
basic set of conditions: N cp = 4; N~ ( = 3; R cp < 0.3 cm; \Z cp \ < 6 cm; number 
of hit wires in the drift chambers corresponds to four tracks. For the events 
surviving those cuts kinematic fit was performed with additional requirements: 
0.1 < E np /2E < 0.3; energy and momentum conservation is held to the accuracy 
higher than 10 MeV; minimal spatial angle between any two particles a m in is 
greater than 9 degrees; total energy deposition of charged particles in the second 
calorimeter layer AEi2 is less than 40 MeV 

In addition, events were required to have some fixed kinematic parameters 
due to presence of intermediate particles in the decay: recoil mass of one of the 
photons must be equal to mass; invariant mass of the other two photons - 
equal to 7r° mass; and invariant mass of this tt° together with two charged pions 
- equal to r\- meson mass. These masses were calculated for analysis, but no 
additional cuts were imposed, 

Detection efficiency under these selection criteria was determined from sim- 
ulation to be equal to e = (0.57 ± 0.07)%. 

Currently the data from three scans i*Y9602 -j- F/9604 were processed, cor- 
responding to a total number of produced <j)- mesons ~ 4.1 • 10 6 . 

Only one event satisfying all intermediate states requirements was found 
(Fig. |l5|). Corresponding branching ratio is equal B(<j> — ► rj'j) = 4 • 10~ 5 . 

The statistics of simulated events of main background processes, available 
by now, is too low to estimate their contributions and only upper limit corre- 
sponding to one detected event could be placed: B{<j> — > 7/7) < 1.7 • 10~ 4 at 
90% confidence level. 



3.4.2 Study of electric dipole decays of light vector mesons 

Electric dipole transitions are widespread in atoms, but in mesons, composed 
of light quarks transitions of the type V — > S7 ( 3 Si — * 3 P a + 7), where V is a 
vector meson (p, u>, <fi) and S is a scalar one (a , f a , . . .) are strongly suppressed 
by a phase space factor because even the lightest scalar mesons are quite heavy 
(M~lGeV). 

Quark structure of scalar mesons is not conclusively determined yet. In a 
standard two-quark description a a f a mesons look like scalar analogs of p- and 
u>- mesons: 

a = and f = or as. 

Better chances has now a four-quark model: 

a = ^^ss and f a = ?^±Mss §§]. 

Yet another model was suggested - kaon molecule, according to which, a a 
f are bound states of two K- mesons (Kj, |ps|). 

And finally, real states can include both ordinary and exotic contributions. 
The probability of transition 4> — > S7 strongly depends on quark structure of 
mesons involved in the decay. Thus the idea emerged |2j| to use radiative decays 
as a probe in order to clarify the structure of a - fo- mesons. The results of 
estimations of branching ratios in different models and the experimental values 
are listed in Table |5[ The most important result is that four-quark model 



11 



predicts branching ratios almost an order of magnitude higher than two-quark 
and kaon molecule models. 



Table 5: Expected values of scalar meson decays branching ratios from the work 
[ pi| p5[ and experimental results. 





qq (ss) 


mm 


KK 


1 


2C 


] 


CMD-2 [ 


21 1 


This work 


B{<t> - fol -» *-<V° 7 ) 
B((f> — > a 7 — » rj-K°^) 


5 • lO" 5 
8 • 10" 6 


10" 4 
10" 4 


10~ 5 


< 2 • 10 _;i 

< 2.5 ■ icr 3 


< 7 ■ icr 4 


(4.7 ± 1.0) ■ 10" 4 
(1.3 ±0.5) ■ 10~ 4 



The first time a search for electric dipole transitions was carried out with 
ND detector |2(|, where upper limits at a level of ~ 10~ 3 were established. At 
present the search for the decays of the kind of 4> — > 57 is being performed at 
VEPP-2M collider by two detectors: SND and CMD-2. A new <jy factory is 
close to commissioning in Frascati where the search for such decays was 
stated as one of the main goals. Another onslaught on neutral decays of <p- 
meson is being prepared at a CEBAF photon beam |3(i[| . So, during next few 
years one can expect, that situation with these decays will be greatly clarified. 

Evidence of the <f> — > 7r°7r°7 decay 
When studying the 

e+e" -> cj> -> 7r°7r°7 (5) 

channel, one should be aware about significant background contributions, 
coming from the following processes: 

e + e~ — > <j> — > rjj — > 37r°7, (6) 
e+e- -> Luir° -> 7r°7r°7, (7) 
e+e" -► ptt° -► 7r°7r°7, (8) 

e + e~ — > (p — > KsKl — » neutral particles. (9) 

Fig. [l6| shows expected recoil photons spectra in the processes (|J), (0), (0). 
One can see, that for £L < 250 MeV contributions from the processes (@, §|) 
are significantly suppressed. The processes (H ^) contribute only due to event 
misidentification, for example, when photons merge together in the calorimeter. 
Events were selected satisfying following criteria: N 1 = 5, 0.8 < E to t/2E < 1.2, 
P to t/2E < 0.12, x% < 25. It was also required, that among all combinations 
of photons existed one, in which two photon pairs have invariant masses within 
the interval m^o ± 30MeV. 

The processes (^|, ||) were suppressed with the help of a special parameter, 
checking if the transverse profiles of photons look like those of isolated ones 
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(xl < 5). Then, the events with a mass of 7r°7 -system lying within the 720 
- 840 MeV interval were discarded in order to reject contribution from the 
process (pp. The ir°ir° invariant mass spectrum of the remaining events is shown 
in Fig. O. Further analysis of the angular distribution revealed, that it agrees 
with S- wave production mechanism of 7r°7r°- system. Detection efficiency for 
the process (^) under described selection criteria decreases linearly with the 
increase of m,.,. from 14% at m w <, w „ = SOOMeV to 6% at 975MeV. This 
dependence was obtained for the process @ simulation with two ir°- mesons 
in S- wave with radiative correction taken into account. The number of events 
and expected background contributions are listed in the Table || 



Table 6: Number of the selected events e + e — > 7r°7r°7 with m n o n o > 800MeV 
and estimated background contributions. 



iV0 


3.9 ■ 10° 


(f> — * 7r°7r°7, experiment 


45 


<fi — > rj-y, simulation 


5 


<j> — > K$Kl, simulation 


< 6 


<fi — » pn°,LOTT° — ♦ 7r°7r°7 , simulation 


1.4 


B((f> -» 7r°7r°7) 


(1.1 ±0.2) ■ 10" 4 



After background subtraction, in accordance with Table ||, 38 ± 7 events 
remain, from which the following estimate for the branching ratio could be 
derived: B{<j> -> 7r°7r°7) = N„o^/N^ ■ e = (1.1 ± 0.2) • 10~ 4 , for 7r°7r°- system 
mass above 800 MeV (only statistical error indicated). So, in this experiment 
we observe the decay (|5j) with S- wave production of tt°- pairs. The comparison 
of 7r°7r° invariant mass spectra, depicted in Fig. with model estimations in 
Fig. ^ shows good agreement with a four-quark model predictions |^5| , 

Fitting, using formulae of |2q] , gives the following f a - meson parameters: 

m fo = (950±8)MeF; 

9\kkI^ = (2.3±0.5)Gey- 2 ; 

9$ omc /4*= (0.4±0.1)Gey- 2 ; 

B{<p / O (980) 7 ) = (4.7 ± 1.0) ■ 10- 4 . 

The latter value was obtained within a framework of the 4-quark model 
(|24 |2q]). Also taken into account was the relation B(f — » 7r + 7r _ ) = 2B(f a — > 
7r°7r°). Systematic error in the branching ratio is estimated to be about 20%. 
To improve it, one should account for the background from the processes ([| ^ 
more precisely (look Table |^) . 

Evidence of the <fr — > ?77r 7 decay 
The main background for the process 

e+e" -> ?77r 7(?7 -> 27, tt° -> 27) -> 57 (10) 
comes from the reactions (|), (|), (|). 

The following event selection criteria were used: N 7 = 5, d min > 27°, 
B 4ot > 0.8 • 2E , X% < 20 ' Xm < 25, x 7 < 20, A^C25 < 6. In the table | the cor- 
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responding detection efficiencies and number of events are listed. Experimental 
statistics here is equal to N$ — 6.5 • 10 6 

Table 7: Detection efficiencies and expected number of events for the main 
processes 





Efficiency 


Number of events 


Experiment 




283 


e + e~ — » rfn-y 


e = (8.7 ±0.3)% 




e + e~ — » 777 


(1.3 ±0.1) ■ 10- J 


109 ±9 


e + e - - Ks K L 


5 ■ 10" b 


10 ±10 


e + e~ — > ujtt 


0.4 ■ 10~* 


85 ±4 



After background subtraction we obtain the total number of events N(cf> — > 
7/7r°7) = 79±21, which corresponds to a branching ratio: B^ir ^) = N/N^-e = 
(1.40 ±0.35) ■ 10" 4 . 

The energy dependence of the visible cross-section was fitted by the sum of 
process ([l0|), background process 4> — > 777 and non resonant contribution from 
the process e + e~ — > um°. The dependence of likelihood function on B(cj) — > 
777r°7) is shown in fig. |l8|. The optimal value is 

B{(j) -> ?77r 7) = (1.3 ± 0.5) ■ 10~ 4 
The mass spectrum of 7y7r°-system is shown in fig. |i~9] 

We also searched for cf> — > rjw 0r y decay in 77 — > 7r + 7r _ 7r° and 77 — > 37r° decay 
modes, but due to higher background and lower efficiencies only upper limits 
for (f> — > 777r°7 branching ratio were obtained: 

B(<j> 777r°7) < 1.5 • 10" 4 (ry -> TT+vr-vr ), 

B(</> ?77r 7) < 3 • 10~ 4 (r? 3tt ) 
at 90% confidence level. 

3.4.3 Search for the <f> — * KsKg-f process 

The 4> ~ * K$Ks"f decay also represents an electric dipole transition, the width 
of which can be expressed as M 2 w 3 , where M is a matrix element and u> is 
a photon energy. The KsKs final state may be produced for example in the 
decays of light scalar mesons a o (980) or / O (980). The expected branching ratio 
of the cf) — * KsKs^i decay, obtained on the basis of the theoretical work |^f| and 
table values on a o (980) / O (980)- mesons |3!| , can lie within a wide range of 
10~ 7 -i- 10 -9 . Currently there are no available experimental data on this decay. 

The search for <fi — > KgKg^f decay was performed in the decay mode <fi — * 
KgKg'y — ► 27r°27r°7 — > 97. The events containing 8 or 9 photons were analyzed 
with the requirement, that invariant masses of four different photon pairs in 
an event must be close to a mass of tt° (m 77 = 110 ± 160 MeV). Then was 
required, that of six possible tt°tt°- pairs at least one has an invariant mass 
close to that of a - meson (430 ± 560 MeV). To suppress background from 
(p — > 777 — * 37r°7 — * 77 decay, the energy of the most energetic photon in an 
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event was required to be less than 360 MeV. In addition the following constraints 
were applied: E tot /2E = 0.8 1.05; P to t/2E < 0.2; x% < 50. 

With such selection criteria the detection efficiency, obtained using simulated 
events of the process under study, varied from e\ = 0.4% up to 0.7%, depending 
on the recoil photon energy oj = m^—mK s K s = 0—24 MeV. The SND threshold 
for photons is about 10 MeV. The detection efficiency with such a threshold is 
equal to (0.7 ±0.1)%. During processing of experimental data from the scans 2, 
3, and 4, total of N a = 16 events were found. All of them could be attributed to 
a background from tfi — > KsKl decay. Thus, the result can be expressed only 
in terms of an upper limit: B{<j> -> K s K sl ) < 2^/W / N^ex < 3.2 • 10~ 4 , at 95% 
confidence level. 



3.5 Rare nonradiative decays of (p- meson 
3.5.1 The process e + e~ — > lutt° — > it + -k~it tt 
Reaction 

e+e~ -> TT + TT-TT°TT° (11) 

was studied in the </>- resonance region earlier [^0| . This process proceeds purely 
via u)ir°- intermediate state and its cross section is close to 10 nb. The main 
contribution comes from p(770)- meson. In the </>- peak region the possibility 
exists to detect contribution from the cj> — ► ujtt° decay, which may reveal itself 
as a narrow interference pattern on a smooth energy dependence of the p(770) 
contribution. Probability of this decay was estimated to be at a level of 5 ■ 10~ 5 

For analysis the events were selected with two charged particles and four 
photons. The energy-momentum conservation and presence of two ir°- mesons 
(Xm < 30) were required. To suppress background from the process e + e~ — > 
K + K~ , restrictions were imposed on specific ionization losses in the drift cham- 
ber. Background from the (p — y m - * vr + 7r _ 7r°7 decay was rejected by a require- 
ment, that the most energetic photon in the event must have an energy less than 
360 MeV. The <fi ~ > tt + tt^tt° decay can also fake the process under study if some 
stray photons are detected. To suppress it, an additional requirement was im- 
posed: if the invariant mass of the two most energetic photons is close to that 
of 7r°- meson, then invariant mass of the other pair must differ from it. 

Event detection efficiency under these selection criteria, estimated using sim- 
ulation, is equal to e = 20%. In the recoil mass spectrum of ir°- mesons, shown 
m Fig. HJ the u>- meson peak is clearly seen. On the other hand, presence of a 
nonresonant contribution is an evidence of some background with intermediate 
states other than ojtt°. To account for this background the special subtraction 
procedure was developed. It is based on background estimation in the kinematic 
regions not containing the process under study. 

The energy dependence of the cross section was approximated according to 
the formula (|^), with Z, the complex interference amplitude, varying, according 
to theory, in the limits: Z = [(11 -=-34) — 24)] • 10~ 2 , and decay probability, 

described by an expression: 
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B{(f> ^ ujtt°) = a (m^)\Z\ 2 /a^. 

For approximation the detection cross section was expressed as in (|l|), and 
the result of the procedure is shown in Fig. ^ll The following values of param- 
eters were obtained: 

a Q = 6.4 ± 0.5 ±1.0 nb - reaction ( pi] ) cross section via wk°- intermediate 
state, 

ReZ = (10 ± 7) • 10" 2 ; ImZ = (0 ± 6) • lO" 2 ; x 2 /N D = 7/9. 

Systematic error of ~ 15% is caused by inexact simulation and background 
subtraction error. As could be seen from Fig. |l]and the values presented above, 
it is possible to establish only the following upper limit of the decay branching 
ratio: 

B((j> -> ujt: ) < 5 • 10~ 5 at 90% confidence level. 

The contribution of non- lutt° intermediate states is a subject of further 
investigation. 



3.5.2 The e+e — > p + n process near <fi- resonance 

Measurements of lepton widths of light vector mesons p, to, tf> at electron- 
positron colliders were usually performed through summation of cross sections 
over all decay channels. Study of the process 

e+e" ->■ fj+fT (12) 

in the vicinity of the c\>- peak gives the possibility of direct measurement of = 
B{4> — > fi + Expected value of B^^, calculated from fj, — e- universality, is 
equal to: 

B^ ~ B ee = (3.00 ± 0.06) • 10- 4 (13) 

The world average value, taken from PDG (1996) p9f : 
= (2.48 ±0.34)- lO" 4 . 

The cross section of the process ([l2]) in the vicinity of the <fi- resonance can be 
expressed (in linear approximation) in the form (^), where Z = 3y/ B^^B^e 113 /a, 
and (3 is an interference phase. 

Events were selected according to following conditions: N cp = 2; A^ 7 = 0; 
A?? < 20°; Aip < 10°; |Z 1)2 | < 8 cm; R < 1 cm; & min > 45°. 

The main sources of background were the following: 

1) e + e~ — > e + e _ , was rejected using e/ir- separation; 

2) e + e~ — > tt + tt~ , was suppressed using the outer anticoincidence system of 
the detector. 

3) cosmic muons. These events are not peaked in time, measured by outer 
scintillation counters, with respect to a beam crossing. On the contrary, the 
events of the process (1) are strongly peaked with a RMS of a ~ Ins. This per- 
mits to reject main part of cosmic muons. Remaining background is suppressed 
by limiting maximum distance from charged particle tracks to the beam axis. 

The detection efficiency for the process (1) is described by the following 
formula: 

= £ sel ' £i3 ' £up • £(i • £flt — 0.32, 
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where e se i — 0.75 is a detection efficiency in the conditions, listed above; 
= 0.55 is an acceptance in polar plane; 

Sup = 0-99 is the efficiency of anticoincidence counters; 

e M = 0.82 probability for a muon to trigger outer system; 

£flt = 0.97 in the second scan, and 0.91 in the third and fourth scans — 
first level trigger efficiency. 

At first, fit was performed with an interference phase as a free parameter. 
The result was j3 = 0.13±0.11. Then the value of f} was set to zero and the fit 
was repeated. The results are shown in Fig. 22 and in Table §. 



Table 8: The results of cross section fitting for the e + e — > fi + \x process 



Experiment 
name 


<r ,nb 








FI9602 
FI9604 
FI9602&FI9603&FI9604 


96.7±1.6 
95.9±1.2 
96.7±0.9 


-0.05±0.21 
0.13±0.12 
0.09±0.09 


2.8±1.4 
2.2±0.9 
2.4±0.8 


9.2/6 
12.2/9 
29/26 



The errors indicated in the Table || are all statistical. Possible sources of 
additional systematic errors are contributions from other decays of <f>- meson 
(for example: <f> — > K + K~, pir, ...); systematic errors in simulation, inadequate 
treatment of radiative corrections. In this preliminary analysis the systematic 
error in <j was estimated to be about 5% and the error in — 20%. Let us 
finally list the main results: 

a a = 96.8 ±0.9 ±5.0 nb, 

B w = ( 2.4 ± 0.8 ± 0.6) ■ 10- 4 , 

= y/B ee • B MM = (2.74 ± 0.44 ± 0.6) • lO" 4 . 

3.6 Rare decays of 77(550)- meson 
3.6.1 Upper limit of the i] — > 2tt° decay 

CP-violation is one of the most intriguing puzzles in the particle physics. Up to 
now the only place, when CP-violation was observed was a neutral kaons system. 
Standard Model (SM) predicts much more pronounced effects in decays of B- 
mesons. 

Flavor conserving CP- violating effects, like 77 — > 27r decay, within a frame- 
work of SM are very weak: B(rj — ► it + it~) < 2 ■ 10~ 27 |Hl p2[ . But some other, 
different from SM, models exist in which such effects are many orders of mag- 
nitude larger. B{r) — > it + it~) < 10~ 15 -j- 10 -16 In any case, search 
for such unusual effects can produce interesting and unexpected results. 

The only existing experimental upper limit — B(rj — > 7r + 7r _ ) < 1.5 • 10 -3 
was reported in the work ]3q| . Recently this result was confirmed by CMD-2 
detector: B(r\ — ► 7r + 7r _ ) < 2 • 10 3 [^TJ] . In this work search for the process 
?y — > 7r°7r°, with ?7-s produced in the reaction e + e _ — * 4> — > rjj was carried out. 
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Events for analysis were selected with two ir°- mesons and a photon, satisfy- 
ing energy-momentum conservation. The estimated detection efficiency is equal 
to 8.6%. In Fig. photon recoil mass spectra in the process e + e~ — > <fi — > 777 — > 
7r°7r°7 are shown. Fig. p3| a depicts an expected spectrum in the process being 
searched for; Figs. |23|b,c show spectra for the main background processes, and 



Fig. 231 — experimental data. Unfortunately the number of simulated events 
is not yet sufficient to determine contribution from <f> — » KgK^ (statistics at 
least equal to the experimental one is needed, i.e. ~ 4 • 10 6 (j>- mesons). As it 
is seen in Fig p3|d, there is no indication of the 77 — > 7r°7r° decay and thus, only 
upper limit can be placed: 

B(i] -> tt°tt°) < 6 • 1(T 4 at 90% confidence level. 

3.7 Nonresonant processes in the <p- resonance region 

3.7.1 The e + e — > 37 reaction 

This process is interesting from the point of view of QED testing near (f>- res- 
onance and as an important source of background in the studies of the decays 
4> —> f]1i 7r °7 ^ 37. 

Events were selected according to following criteria: N 7 = 3; i) m in > 27°; 
0.75 < E tot /2E < 1.1; P tot /2E < 0.1; X * < 15; x| < 30. 

The Dalitz plot for the selected events is shown in Fig. ^4|. Regions of 
4>- meson decays (f> — > 777, 7r°7 — > 37, are marked and events from them were 
excluded from further analysis. The energy dependence of the visible cross 
section for the selected events is shown in Fig. ^5[ The estimated cross section 
at 510 MeV is equal to <jmc — 2-03 ± 0.07 nb, while the experimental value is 
20% smaller: a exprm = 1.56 ± 0.16 ± 0.20 nb. The cause of such discrepancy 
is under study. In Fig. the spectrum of the least energetic photon in an 
event is presented. One can see, that the largest difference between theory and 
experiment lies in the region E m i n /E ~ 0.5. 

The e + e~ — > 37 process is a source of background for neutral decays of 
vector mesons and, on the other hand, its cross section is precisely known. Its 
investigation is important for understanding of systematic errors in other neutral 
processes and accuracy of luminosity estimations. 

3.7.2 The process e + e — > e + e 7 and search for 7/ — > e + e decay 

The process e + e~ — > e + e~7 is important for us due to several reasons: 

1. this process can be used for testing of QED and should be taken into 
account as a correction for Bhabha scattering, in order to achieve accuracy 
of luminosity measurements of < 1%. 

2. this process is a background to different hadronic processes, i.g. — > 
7r°7, 777 — > e + e - 7. 
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This process was studied earlier in several experiments: with OLYA detector 
IH with an integrated luminosity of AL = 1.5 pb' 1 , at ADONE collider |$7| 
with AL — 0.12 pb^ 1 . 

The results of this work are based on data corresponding to integrated lu- 
minosity of 1.8 pb^ 1 . During data processing the subset of events with Dalitz 
pairs was also included. In such events the invariant mass of e + e~- pair is very 
small and the whole process looks like a two-photon annihilation in which one 
of photons is slightly out of a mass shell, and it subsequently decays into e + e~- 
pair. 

The e + e~7 final state is favorable for the search of the decay e + e~ —>(/)—> 
ryy — > e + e~7. The 77 — > e + e~ decay is strongly suppressed in comparison with 
■q — ► 77 by a hclicity factor (rn e /m v ) 2 and additional QED suppression factor 
a 2 . As a result the expected branching ratio is about 5 • 10~ 9 |f38[ . The similar 
decay tt° — > e + e~ was observed at a level of ~ 10~ 7 . Current upper limit of the 
77 -> e+e" is equal to 2 • 10~ 4 @ 

Event selection was done using the following main criteria: N cp — 2; — 
1,2; d mm > 36°; R < 0.5cm; \Z\ < 10cm; E tot /2E > 0.8; P tot /2E < 0.15; 
A(p ee > 5°; Xe < 40; Af2 ee < 50° — spatial angle between particles for Dalitz 
pairs. 

Fig. ^ shows the photon energy spectrum and in Fig. ^8| the distribution 
over angle between electrons is shown. There is also a region indicated, where 
the angle between electrons is very small (Dalitz pairs). The number of events 
in this region is in agreement with an expected one. 

Energy dependence of the e + e~ — > e + e~7 cross section was approximated 
using the following expression: 

O-vis(E) = Ei(Jo(%-) 2 + SO-31T- 

The first term is responsible for the process under study, while the second 
describes background from decays of 0- meson, especially <j) — » 37r. It was esti- 
mated from background simulation, that e = (0.10 ± 0.03)%. The £1 coefficient 
was determined as a ratio between the number of simulated events, satisfying 
selection criteria and the total number of simulated events. The simulation was 
conducted with the following constraints: i? m j„ > 27°, E 1 > 5 MeV, and corre- 
sponding total cross section under such constraints is equal to <j Q = 570 ± 15 nb, 
and e\ was found to be equal to (3.5 ± 0.1)%; 

The following results were obtained for approximation parameters: a Q = 
550 ± 10 nb; e = (0.08 ± 0.05)%; X 2 /N D = 10/11. Only statistical errors are 
indicated. Additional systematic error is estimated as 10%. 

The 77 — > e + e~ decay could reveal itself as a peak in the photon spectrum in 
Fig. [27] at E 1 = 360 MeV. The detection efficiency, obtained from simulation 
is £2 = (40 ±1)%. After fitting of the spectrum the following upper limit was 
obtained: B(i] — > e + e _ ) < 9 ■ 10~ 4 . Although this limit is higher than that of 
PDG Table, increase in statistics and taking into account the energy dependence 
of the spectrum shape will permit us to improve this result. 
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3.7.3 e+er -> lutt" -> 7r°7r°7 

There are several reasons of interest to the e + e~ — > W7r° process. First, in 
this process may manifest themselves radial excitations of p- meson. Second, 
this process contributes into the total cross section of e + e~- annihilation into 
hadrons. Third, its cross section is connected with the branching ratio of r- 
lcpton decay into wkv t . 

Within the Vector Dominance Model (VDM) cross section of this process is 
determined by g pU i-n constant, which also controls other processes: to — > 37r, 7r°7; 
p — > 7r°7; 7r° — * 27, etc. 

In our case the possibility arises to measure this constant at y/s = 1 GeV, 
which can help to determine boundaries of applicability of simple VDM model. 

The process e + e~ — > ujtt in the energy region 2E > 1 GeV was studied 
earlier in the experiments with ND detector at VEPP-2M in two main decay 
modes of u>- meson: u> — » tt°j J42| and ui — ► tt + tt~tt° |2^, 

Some indirect results on this reaction were extracted from r- lcpton decay 
spectra 

In this paper data on the process e + e _ — > oj7r° — ► 7r°7r°7 in the energy 
range 2E 1 = 1000-^ 1034 MeV are presented. Results are based on experimental 
statistics from FI96 SND experiment corresponding to integrated luminosity of 
1.9 pb^ 1 or 4 million <fi- mesons produced. 

The main event selection criteria are: 7V 7 = 5, E to t/2E = 0.75 -7- 1.25, 
d min > 27°, P tot /2E < 0.1. 

Then additional requirements were imposed: xti < 40; x 7 < 20; two photon 
pairs with ± 30 MeV which could be produced by intermediate tt°- 

mesons must be found in the event; 720 < m^o^ < 840 MeV. 

Selected events originate mainly from the process under study: e + e~ — * 
lot:" — > 7r°7r°7. However they contain some background from the decay (f> — > 
777 — > 37r°7. The detection efficiency was estimated by simulation. For the 
process under study it is equal to e\ — 21% and for the background — £2 = 
(1.8 ± 0.1) ■ 10~ 3 . As for another background source: (f> — > KsKl — > neutral 
particles, the available simulation statistics is not sufficient to estimate its con- 
tribution (minimum 10 6 events are required, while only 6 • 10 4 are currently 
available). 

In Fig. ^ the distributions over invariant mass of 7r°7- system in experimental 
and simulated events are shown. It is clearly seen that most of events contain 
intermediate u>- meson. In Fig. ^ of two possible 7r°7- pairs in an event only 
one with the mass closest to that of oj- meson is histogrammed. The number of 
experimental events in the u- peak region is equal to 304, while expected back- 
ground from the decay </> — > 777 is 92 events. The background from 4> — > KgK^, 
according to available simulation statistics, is not more than 20 events. This 
background could be further suppressed by additional constraint on invariant 
mass of ir°ir°- system to exclude the K3 — > 2ir° decay region. 

The energy dependence of the e + e~ — > lot: cross section was approximated 
using parameterization (||) with B(ui — > 7r°7) = 8.5 ± 0.5%, and taking into 
account contributions of background processes — > KsKl,t]j into crs(-Bi). In 
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this approximation due to small statistics the contribution from tf> — > pn, ujt: — > 
ir°ir°j was neglected and thus the energy dependence of the cross section was 
described by a simple linear function. As a result of fitting of the visible cross 
section (Fig. ) the following results were obtained: 

o (m^,) = 7.5 ± 1.5 ± 0.5 nb — the cross section of e + e~ — > ujir° — > n°TT°j 
in the 0- meson peak (previous ND'1991 result (p0[)was cr o (m0) = 8.7 ± 1.0 ± 
0.7 nb); o'(m$) = 9 ± 6 MeV~ x — relative slope of the cross section at the 
same energy; P(x 2 ) = 12.3/9. 

The <7 pw7r value in an assumption of p-meson dominance is equal to: 

9punr (experimental) — 18.1 ± 1.7 GeV^ 1 . 

in VDM this constant equals to: 

9punr (estimated) = 14.3 GeV^ 1 . 

In the ND experiment (2(| [|2| the cross section was larger, but it does not 
contradict the new SND result. The SND data confirm conclusion of the ND 
experiment that pure p(770) production cannot describe the data and to take 
into account production of higher resonances p', a/, <f>' is necessary. 

3.7.4 Production of A- barion 

When electron beam travels through the collider beam pipe the following process 
may occur on nuclei of residual gas atoms: e ± A — > e ± A — > e^Nit pEfl . The 
A(1232) is a barion state with quantum numbers: I(J P ) — |(| + )- At a beam 
energy of E = 510 MeV the total cross section of this process ~ 3ub. A(1232) 
decays into the following final states: A + — > p7r°, nir + A° — * pir~ , mr°. 

An attempt was taken to observe this process with SND detector and to 
study it as a possible source of background at future cj)- factories. 

The process was studied in the: 

e ± p -> e ± A+ -> e ± P Tr° (14) 

channel because it has a distinct signature: two photons from 7r°- meson de- 
cay, slow proton track with high dE/dx in the drift chamber, another track 
of e^, zero total transverse momentum and large longitudinal momentum - 
P = 510 MeV/c. 

The events with N cp = 2, N 7 = 2, E tot /2E < 0.5, \Z\ < 10 cm, etc. were 
selected. Then kinematic fit was performed. 

As a result of experimental events processing in the energy region of cf>- peak 
with integrated luminosity of AL — 0.5 pb^ 1 , 36 events were found. Photons in 
these events have an invariant mass coinciding with that of tt°- meson within 
±15 MeV. These events are uniformly distributed in Z- coordinate within 
±10 cm from the interaction point and show no peaking there. Among two 
charged particles, one with lower specific ionization was referred to as an electron 
and another - as a proton. In Fig. [3(j]the specific ionization in the drift chamber 
as a function of particle momentum is shown for particles after kinematic fit for 
experimental and simulated events. One can see, that protons are well separated 
from electrons. The A- barion mass spectrum in Fig. |3l] shows good agreement 
between experiment and simulation. 
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The background from <j>- meson decays is still questionable because of insuf- 
ficient simulation statistics, but it is obviously not high. It can be seen from 
absence of peaking of the events near the beam interaction point. 

Detection efficiency was estimated from the simulation and was found to be 
e = 2%. Given the pressure of residual gas of about 3 • 10~ 9 Torr, mean beam 
current of 15 mA, duration of the experiment 7 • 10 5 s, and taking into account 
residual gas composition of 50% - H 2 , 30% - CO, 20% - C0 2 f4§, one should 
expect N ~ 4000 events of A- barion electroproduction at 20cm interaction 
region length. Taking into account the detection efficiency, the number of ob- 
served events is expected to be equal to 30 (|l4|). It agrees well with N Q = 36 
events detected in experiment. We should mention, that there is a considerable 
uncertainty in the pressure and chemical composition of residual gas. Presence 
of CO2, N2, etc. may lead to a considerable increase in events number. 

In conclusion, we can say, that event rate of this process ~ 6 • 10~ 3 is too 
small to affect other processes under study, but at future (jy- factories it should 
be taken into account (|lj) as a possible source of background. 
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4 Recent results 



In August 1997 just before printing of this preprint the following new upper 
limits for processes with multiphoton final states were obtained: 

• S(0^2tt°) <4-l(T 5 

• B{4> -> ?77r ) < 5 ■ 1CT 5 

• B((j) -+3tt°) <0.8-10~ 5 

The event selection criteria were similar to those applied in the studies of 
the processes described above, like <f> — > r/j — > 37r°7 or <fi — > rjir j. 

Also was performed search for K$ — > 37r° decay. The event selection criteria 
were adjusted for 4> — * KsKl events with i^g — > 37r° decay and with completely 
undetected K^. No indications of Kg — ► 37r° decay were found and the upper 
limit B(K S -> 3tt°) < 1.1 • 10~ 4 was placed. 



5 Conclusions 

The SND detector continues data taking at VEPP-2M collider. In 1997 the 
MHAD97 experiment in the energy range 2E from 980 up to 1380 MeV and 
total integrated luminosity of 6.3 was carried out. The work is under way 
on improvement of energy and spatial resolution of the detector. About 50% of 
experimental statistics of the FI96 experiment, collected in 1996 in the energy 
range 2E = 984 1040 MeV were analyzed. New preliminary results on <fi- rj- 
mesons decays were obtained (Table ^) . 

Cross sections of the following nonresonant processes were measured in this 
energy region: 

a (e + e- -> ujtt°{lo -> tt+tt-tt )) = 6.4 ± 0.5 ± 1.0 nb; 

a (e+e- -> lutt°(uj -> 7r 7 )) = 7.5 ± 1.5 ± 0.5 nb; 

a (e+e~ -> = 96.8 ± 0.9 ± 5.0 nb; 

er (e+e- -> 777) = 1.56 ± 0.16 ± 0.20 nb, at d mm > 27°; 

a a (e + e- -> e + e~i) = 550 ± 10 ± 55n6, at d mm > 27° £ 7 > 5 MeV. 

Electroproduction of A- barion on residual gas nuclei was seen: e^A — ► 
e ± A -> e ± Nir. 
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Table 9: SND preliminary results on 0(1020), 77(550) , Ks decays, compared 
wi th CMD-2 g|, H results and data from PPG @ 



Decay mode 


OMR 




ir L)vj(iyyb ) 


<j> — > ?77 


/-1 nri 1 n _l_ n n*7\07 
(l.oU ± U.Ub ± U.U/J/o 


(1.2o±U.lU)/o 


(1.2b±U.Ub)/o 


17 — * 77 


^1.0/ it U.U4 it U.UoJ/o 






// — ► OvT 


9^ -t- n n 1 } 4- n 

l^l.ZO It U.UO It U.UO ) /o 






— — ^ 


^1.1)^1 It U.UO It U.UO j /o 






—* 7T 7 


^u.±o n u.uuu it u.uuo^/o 




in 1 qi 4_n m ^"l^ 


ip — > r; ^aoo;7 


< 1 7- 10~ 4 


< 2-10 - 4 


< 4.1 ■ 10 — 4 


> KsKs'y 


< 3.2-10 - 4 






— + 7r <) 7T°7 


H 1 4_n 9V10 - 4 
^ j- . ± in ' j.u 




< 1 ■ 10 — 3 


— ► 7T^7r — 7 




< 1.510" 5 


< 7-icr 3 


-» M+/i"7 


— 


(2.3±1.0)-10" 5 


— 


— » r)7r°7 


(1.3 ± 0.5)-10~ 4 


— 


< 2.5-1CT 3 


-» /07 


(4.7±1.0)-lCr 4 


< 1 - 7-10" 4 


< 2-icr 3 


-> p 7 


— 


< 310" 4 


< 7-icr 2 


-> tfg^ 


0.31±0.02 


0.325±0.010 


0.341±0.005 


K s -> 3tT° 


< 1.1-10 -4 


— 


< 3.7-10" 5 


-> ft+fl - 


(2.4±1.0)-l(r 4 


— 


(2.48±0.34)-lCT 4 


— » e+e — 




(2.84±0.0b)-lCr 4 


(3.00±0.06)-10" 4 


— > r;e + e — 




(i.i±o.5)-icr 4 


(i.3±°; 8 6 )-i(r 4 


— > Tr^e^e - 




(l.l±0.8)-lCr 5 


< 1.2-10 -4 


— ► o;7r 
(7r+7r-7r°7r°) 


< 7-1CT 5 






0— OT^7T — 7T^7T — 




< 3-1CT 5 


< 8.710" 4 


-> 2tt° 


< 4-10" 5 






-» 3tt° 


< 0.8-10" 5 






— » r)7T° 


< 7-1CT 5 






77 — * e + e — 


< 9-10" 4 




< 3-10" 4 


7) — > 7r 7r° 


< 7-l(T 4 






77 — > 7r~*~7r~ 




< 9-icr 4 


< 1.5-10" 3 
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Figure 1: Average luminosity (integrated luminosity at each energy point di- 
vided by corresponding live time of data taking) in the experiments with SND 
detector as a function of the beam energy. 
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Figure 2: SND total integrated luminosity accumulation schedule in the period 
from 01/13/1997 to 06/13/1997. 
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Figure 3: Distribution of e + e~ — > 77 events over normalized energy deposition 
in the calorimeter. Points — experimental data, a exp . — 5.0%; histogram — 
simulation, om.c. — 4.2%. 
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Figure 4: Distribution of the events over 7r°7 invariant mass in the e + e 



lutt° 7r°7r°7 process; a 



experiment, b — simulation. 
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Figure 5: Acollincarity angle distribution in azimuth plane in the e+e — > e+e 
events. 
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Figure 6: Polar acollinearity angle distribution in the e+e — > e+e events. 
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Figure 7: Distribution over the distance from a track to the beam axis AR in 
the e + e~ — > e + e~ events. 
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Figure 8: Ratio of the luminosities L 77 /L ee in the FI9602 scan 
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Figure 10: Cross section of the e + e — > ip — > KgK^ process. 



33 



^2000 



1500 



1000 



500 

















a =16.1 
1 a =74.9 

mc 






















1 





450 



500 



550 



600 



Figure 11: r\- meson invariant mass in the e + e — ► <p — > 777 process. 
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Figure 12: Cross section of the e + e — > ^ — * 777 process. 
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Figure 13: Recoil mass distribution for the most energetic photon. 
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Figure 14: Distribution over the reconstructed mass of the r\- meson; points 
with error bars — experiment, a exp — 16.6 MeV; histogram — simulation, 
a sim = 16.9 MeV. 
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Figure 15: m rec . 7 vs. 
7r + 7T~7r + 7r~37 decays. Circles 
selection criteria. 
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Figure 16: Recoil photon spectrum in the process e+e — > 7T°7r°7. 
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Figure 19: Measured mass spectrum of rjn system. 
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Figure 20: Recoil mass spectrum of ir°- mesons in the e + e 
process. Top — simulation; bottom — experiment. 
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Figure 21: Detection cross section and fit parameters for the e + e 



7T + 7r 7r°7r° process. 
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Figure 22: 



cross section. 
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Figure 23: Photon recoil mass spectrum in the e + e~ — ► cp — > 7r°7r°7 process; 
a — X] — > 27r° simulation; b — W7r° (cj — ► 7r°7) simulation; c — 77 — > 37r° 
simulation; d — SND experiment. 
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Figure 24: Dalitz plot for the events of the e + e — > 37. process 
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Figure 25: e+e -^37 detection cross section. 
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Figure 26: Lowest energy photon spectrum in e + e — > 37 events. 



n 



n 


- experiment 

- MC simulation 




- E y for Br{y\-^tee) 



3 

1500 
1000 V 
500 




Figure 27: Photon spectrum in the e + e — > e + e 7 process. Peak at i? 7 = 0.7.E 
comes from the possible 7/ — > e + e~7 decay with a branching ratio close to 1%. 
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Figure 28: Distribution over the spatial angle between electrons in the e + e 
e + e~7 process. 
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Figure 30: Dependence of dE/dx in the drift chamber on particle momentum. 
Experimental and simulated data for electrons and protons from the process 
e^p — ► e ± A + — ► e^~p7r° are shown. 
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Figure 31: Distribution over invariant mass of proton and pion in the e ± p 
e A + — > e ± pir° process 
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